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ABSTRACT 


The  exact  formulation  of  displacements  and  resonant  frequencies  of  an  elastic 

w M 

solid  cylinder  has  been  attempted.  Progress  hft****«n  made,  hut  the  principal 
objective  m  not  tea*  reached.  This  report  traces  the  course  of  work  through 
the  early  demonstration  of  invalidity  of  a  formulation  that  was  postulated 
prior  to  the  contract,  and  through  several  assaults  upon  the  problem.  Novel 
solutions  of  the  wave  equation  for  radially  symmetric  displacements  hrra 
discovered.  A  section  is  devoted  to  the  recommended  direction  of  future  work. 
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Section  I 
INTRODUCTION 

"  ,  t.  v 

This  is  the  ^Technical ^Rapert  issued  upon  completion  of  Contract 
AF  1*9(638)  -  11U8.  It  covers  the  period^jjjg  Mgfch  tiMWCgh  lli  December 
1962*  This  contract  with  Vitro  Laboratories  has  been  entitled  "Elastic 
Solids  Vibration  Studies". 

It  has  been  the  objective  of  work  under  this  contract  to  advance  an 
important  area  of  the  theory  of  elasticity.  This  area  is  the  mathematical 
description  of  the  elastic  motions  and  resonant  frequencies  of  isotropic 
solids  that  are  bounded  by  traction-free  surfaces.  Excepting  the  sphere,1 
we  may  say  that  exact  solutions  for  solids  of  arbitrary  dimensions  with 
these  surface  conditions  are  not  yet  known}  other  works  have  either  pre¬ 
sented  (1)  exact  solutions  for  solids  that  are  semi-infinite  in  one 
dimension  or  (2)  approximate  solutions  for  solids  of  arbitrary  dimen¬ 
sions. 


1A.  E.  H.  Love,  "A  Treatise  on  the  Mathematical  Theory  of  Elasticity", 
Uth  edition,  Dover  Publications,  New  York,  19U1*,  pp  278-287* 


Section  II 

THE  ISOTROPIC  SOLID  CYLINDER 

RESTRICTIONS 

The  right  circular  cylinder  of  arbitrary  dimensions  is  the  model  studied 
during  this  contract.  To  reduce  the  mathematical  labor,  While  yet  demonstrat¬ 
ing  the  techniques  of  solution,  attention  has  been  restricted  to  modes  in  Which 
motion  is  symmetrical  about  the  axis.  This  casts  the  problem  in  two  dependent 
variables,  one  more  than  the  sphere  presents  in  its  maximum  symmetry,  yet  one 
less  than  presented  by  the  rectangular  parallelepiped  of  arbitrary  dimensions. 
Quite  arbitrarily  the  work  has  been  restricted  also  by  considering  only  those 
vibrations  Which  are  symetrical  about  the  midplane  of  the  cylinder. 

A  BOUNDARY  VALUE  PROBLEM 

The  starting  point  in  this  contract  was  a  set  of  displacement  and  fre¬ 
quency  equations  previously  developed  by  the  writer.  Within  the  scope  of  this 
contract, 'it  was  the  writer's  initial  objective  to  confirm  these  by  various 
tests  before  proceeding  into  a  computational  phase  in  Which  fundamentals  and 
overtones  would  be  found.  This  was  to  be  followed  by  an  experimental  program 
to  verify  frequencies.  However,  as  will  be  seen  in  subsequent  paragraphs,  the 
course  of  the  study  was  sharply  altered  when  the  confirmation  of  initial  e- 
quations  failed. 

The  displacement  and  frequency  equations  resulted  from  a  derivation  which 
formulates  as  a  boundary  value  problem  the  elastic  motions  and  dynamic  stresses 


in  the  isotropic  cylinder.  The  entire  derivation  is  presented  in  Appendix  A, 
and  the  displacements,  resonant  frequency,  and  components  of  resonant  frequency 
are  reproduced  here  as  Eqs.  (l)-(5). 


Radial  and  axial  displacements  are,  respectively 


Ur.  =i[A)}  Pi  ^i(e>A-) 

J  A 


Cos  m;z  + 


A)j+t  hy+i  J,  cos 


-  L  IaJI 


(1) 


(2) 


con- 


Ut  =  <EjiAj)  mj  sin  rr\j  z  - 

Ajj+i  fj  J»  (fj*)  sin.  tnf+,z  ]&'Cu,t 

Where  and  are  arbitrary  constants}  Pj,m^,nij+^  are  propagation 

stantsj  and  on  is  angular  frequency. 

Note  that  the  displacements  are  symmetrical  about  the  axis  of  the  cylinder, 
and  are  also  symmetrical  about  the  midplane  of  the  cylinder.  The  latter  property 
could  as  well  have  been  made  antisymmetric,  in  Which  case  the  displacements  would 
describe  the  flexural  modes. 

Angular  frequency  »  in  Eqs,  (1)  and  (2)  is  related  to  propagation  constants, 
density  p  and  Lame  constants  X  and  by 

-fgr^  =  (ft  -«<;<*>  (3) 

Propagation  constants  are  solved  from  the  following  pair  of  simultaneous, 
transdendental  equations 


ViS  F ( PH*)  f  (pH  -  f  F  ( PH)  »  a  ejx  (  pf_,  +  mj2)  (U) 

0  (5) 


/' 


Where 

and 

THE  STAIRCASE 


•Ui> 

f(q)  -  q  tan  q 


(6) 

(7) 


The  meaning  of  Eq.  (3)  is  illustrated  By  Fig.  1  in  Which  the  two  lines 
inclined  U5*  to  the  p2  and  m2  axes  represent  the  locus  of  points  of  the  same, 
constant  .  The  inner  of  the  two  lines  is  associated  with  the  dilatational 
components  of  Eqs.  (1)  and  (2),  and  these  are  the  terms  in  Aj^j  the  outer  line 
is  associated  with  shear  components.  For  Poisson's  ratio  of  0.3,  the  dila¬ 
tation  and  shear  intercepts  on  either  of  the  axes  stand  in  ratio  of  2  to  7 • 

The  staircase  effect  displayed  by  Fig.  1  is  formed  by  connecting  all 
points  which  are  predicted  by  Eq.  (3)  after  selecting  a  specific,  allowed 
pair  of  Pj,fflfc,  as  for  example  based  upon  the  point  (p02,»^2)* 


-U- 
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It  is  shown  clearly  by  Fig.  1  that  predicted  pairs  of  p  and  *  will 
consist  of  one  positive  and  one  negative  value,  or  both  will  be  positive* 

It  is  Impossible  for  both  members  of  the  pair  to  be  negative.  Since  the 
type  solution  for  displacement  is  a  Bessel-trigonometric  product,  a  negative 
value  makes  one  member  of  the  product  a  hyperbolic  function,  recognisable  as 
a  Rayleigh- type  wave* 

A  MATHEMATICAL  DIFFICULT! 

The  illustrative  point  (p '02#j^2)  is  readily  found  from  Eqs.  (U)  and  (5), 
The  writer  assumed  pc2  -  m©^  •  1  and  found  pQa  »  1.36$,  x^h/2  -  1.0U18,  and 
hence  h/2a  -  0.761*.  By  Eq.  (3),  all  values  of  Pj2  and  n^2  for  j  /  0  were  now 
known.  However,  it  was  at  once  found  that  these  other  values  did  not  simulta¬ 
neously  satisfy  Eqs.  (U)  and  (5)  when  the  latter  were  written  with  J  /  0. 
Therefore  Eqs.  (l)-(5)  do  not  constitute  a  valid  aolutionj  there  are  too  many 
equations  for  the  number  of  unknowns. 


COMPLEX  PROPAGATION  CONSTANTS 

A  coup lex  propagation  constant  is  admissible  in  Eq.  (3)  provided  that 
frequency  remains  real.  If  pj  »  ♦  if ,  and  *  ivj  in  Eqs.  (l)-(5)* 

the  number  of  unknowns  in  the  entire  set  of  equations  is  doubled,  and  the 
number  of  equations  is  doubled.  The  mathematical  impasse  found  with  Eqs. 

(l)-(5)  appeared  to  be  unresolved  by  substitution  of  complex  propagation 
constants. 

A  THIRD  EXPRESSION  FOR  FREQUENCY 

An  alternative  to  the  creation  of  more  unknowns  with  Which  to  validate 
Eqs.  (3)— (6)  is  to  terminate  in  a  few  terms  the  infinite  series  represented 
by  Eqs.  (1)  and  (2).  It  occurred  to  the  writer  that  a  third  wave  equation 
might  lurk  in  the  vector  wave  equation  of  elasticity  in  isotropic  media, 

+  VV-U  -  VX  VX  U  a  (8) 

Theorists  in  elasticity  write  of  one  congressional  and  two  shear  waves  in 
any  general  discussion.  The  writer  had  only  two. 

An  effort  was  mounted  to  see  if  Eq.  (8)  could  possibly  be  converted  to 

=  °  <’> 

or  into  any  3-bracket  operator.  (The  symbol  $  in  Eq.  (9)  is  the  vector 
Laplaclan  operator.)  The  result  of  this  work  was  the  following  expression 

[(  A+2/*;  2*  -  /» ]\m  &  ](>*  “7°  0  (10) 

which  produced  a  third  line  of  constant  frequency  in  Fig.  1.  This  third 
line  made  it  possible  to  replace  the  infinite  staircase  with  a  closed  figure 
involving  six  pairs  of  p  %m  values.  The  series  would  then  be  terminated  in 
6  terms. 
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However)  the  third  bracket  would  not  satisfy  the  equation  of  motion  from 

2  V 

which  it  had  come  —  unless  •  ■  0.  Therefore  the  conversion  of  Sq.  (9)  into 
Bq.  (10)  had  introduced  multiplication  by  zero  and  was  invalid.  This  was  later 
shown  to  be  the  case.  The  direct  conversion  of  Sq.  (9)  into  wave  equation 
factors  can  be  accomplished  by  operating  upon  both  sides  of  Bq.  (9)  with  the 
factor 

Q MVV'  -  (k  +  2.yu)V*V  *  ~/°&] 


whereupon  one  arrives  directly  at  the  usual 

<u) 

In  performing  this  operation  one  needs  the  vector  identities 

VX  vx  vx  fxA  =  -  v ?•  VV-A  (i2) 

and 

W’VV'A  =  YV-&A  =&  W-A  (13) 

It  has  since  been  concluded  that  the  second  shear  wave  is  missing  from 
the  present  cylinder  problem  because  of  symmetry  considerations.  It  would 
take  the  form  of  displacement  in  the  •  direction  both  for  a.  wave  propagating 
radially  and  for  one  traveling  parallel  to  the  amis.  However)  a*  is  totally 
decoupled  from  %  and  in  the  present  simplified  problem. 

ADDITIONAL  SOLUTIONS  TO  WAVS  EQUATION 

It  is  declared  in  Appendix  A  that  the  admissible  functions  for  radial 
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and  axial  displacements  are 

=  A  Cos  m  =  (Jn  tTiuii  (1U) 

Uz  -  S  J*(f*-)$in  n*  e.~LlJ^  =  Uz  &  (1$) 


During  this  contract  discovery  has  been  made  of  other  admissible  solutions 

of  the  wave  equation.  The  simplest  of  these,  with  time  suppressed, 

are 

C  [toA  J0(pA)  L*s  mz  +  f7  Jj  (fA.)  Sih.  Wi  j  (16) 

If  r  +  />*  J«^a)  (17) 

and 

^  --  ■/-(wa)4]  J((fA)  caj  /mi  -f  ^»a  J0(j>A)(cos  mi  -  2wz  ;/»i  mii)| 

(18) 

(/*»  F  ([  (f  x)* +  (**)'■]  Jt(r)  Sin  »z  -  A  J,(fn.)l(^i-Zr^)sin  <*z 

+  2f>xmz  cosmzj 

(19) 


Substitution  of  Eqs.  (l6)-(19)  in  the  equation  of  motion  is  illustrated  in 
Appendix  B. 

It  may  be  seen  from  Eqs.  (lU)-(19)  that  the  general,  two-displacement 
solution  of  the  wave  equation  is 

U,  =  £  £  Aimz2t{[A»«,xj  J  ,((>«■)  +  A  J,Q*)]cosrhz 

^  ■+  ftJ/fm) ]z  sin  m-z  j 


*  o 


sm  mz. 


1  f  *  1 

a  JiCf**-)]?  C W  ftizj 


(21) 


where  Eqs.  (1U)-(19)  are  subsets  of  the  general  solution, 


* 

< 


} 


1 


The  arbitrary  constants  A:2nf2q,****®2n4'l  2q+l  w^’  136  determined  toy 
the  constraints  of  the  problem.  From  this  point  of  view  there  are  seven 
constraints,  namely,  the  two  vector  wave  equations,  the  vector  equation 
of  motion,  and  the  four  zero-stress  conditions  upon  the  flat  and  curved 
surfaces. 

Considerable  time  during  the  contract  period  has  been  spent  without 
significant  progress  in  an  effort  to  establish  these  constants.  One  of 
the  results  is  that  the  constants  Aqq  and  B00  are  independent  of  any  other 
constants,  Which  verifies  the  initial  use  of  zero  order  terms  in  Eqs.  (ll±) 
and  (15>) .  An  insight  which  has  been  developed  during  this  effort  is  that 
certain  simplifications  of  the  matrix  of  constants  may  be  in  order.  Ex¬ 
amples  of  these  are  rejection  of  constants  with  even-even  subscripts,  i.e., 
A2n,2q  *  B2n  2q  "  constoants  with  odd-odd  subscripts,  i.e,,  A2n+i}2q+l  "  - 

B2n+l,2q+l  *  °J  or  of  constants  with  subscript  pairs  in  Which  the  two  indices 

differ  by  more  than  four  or  perhaps  three,  e.g.,  A2n  2q  “  •  *  A2n+l,2q+l  "  0 

2  2 

for  (2n  -  2q)  >  16  or  (2n  -  2q)  >9,  etc.  No  reliable  guidelines  to  such 
rejection  have  as  yet  been  discovered. 

/ 

THE  CALCULATION  OF  PAIRS  OF  (p02,B^,2) 

In  Appendix  B  are  shown  the  three  simplest  forms  of  radial  and  axial 
displacement  (through  increasing  n  and  q)  after  substitution  in  the  equation 
of  motion.  Both  dilatation  and  shear  forms  are  presented.  When  these  are 
subjected  to  the  stress  conditions,  there  results  from  the  pair  containing 
the  arbitrary  constant  A  the  two  transcendental  Bqs.  (U)  and  (5)  as  seen 
earlier.  With  some  surprise  it  was  found  that  Eq.  (5)  again  emerged  from  - 


u 
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the  flat- surface  conditions  when  they  are  applied  to  the  pair  in  arbitrary 
constant  B,  Still  a  third  time  does  Bq.  (5)  evolve  from  the  flat-surface 
conditions  when  they  are  applied  to  the  pair  in  arbitrary  constant  C. 

Although  the  curved-surface  conditions  have  not  up  to  this  time  been 
applied  to  the  second  and  third  simplest  forms  of  displacement,  it  was  felt 
that  the  unfailing  appearance  of  Eq.  (£)  might  have  significance,  that  it 
might  be  saying  that  this  pair  of  equations,  Bq.  (U)  and  (5)  when  restricted 
to  j  -  0,  are  correct  in  predicting  frequencies}  and  that  there  remains  for 
the  future  only  the  hurdle  of  decoupling  from  the  3-0  case,  through  better 
choice  of  functions  for  displacement,  the  pairs  of  transcendental  equations 
for  3/0. 

Accordingly  an  exploratory  computation  of  the  roots  of  Eq.  (U)  and  (£) 
for  3-0  has' been  made.  The  results  of  this  effort  are  presented  tabular ly 
in  Table  I  and  graphically  in  Figs.  2  and  3.  Plotted  points  are  calculated 
values  from  Which  curves  have  been  generated. 

The  curves  in  Fig.  2  are  smooth  whereas  those  of  Fig.  3  suggest  that 
roots  associated  with  more  than  one  mode  have  been  confused.  More  computation 
is  desirable  in  order  to  provide  (l)  an  unscrambling  of  roots  in  Fig.  3}  (2) 
a  wider  range  of  (»o/Po)^<  0  in  Fig.  3}  (3)  and  higher  order  roots  in  both 
Figs.  2  and  3* 


.11. 


TABLE  I 


CALCULATED  ROOTS  OF  EQS.  (3)-(5)  FOR  j  -  0 


”o 

i^h 

A 

Poa 

Po 

2 

31.62 

1.568 

.0761* 

10.00 

1.5U2 

.2383 

U.U72 

1.1*81 

.501* 

3.162 

1.1*23 

.673 

2.00 

1.313 

.936 

l.i*H* 

1.192 

1.152 

1.00 

1.0U2 

1.365 

0.775 

0.9201* 

1.509 

0.707 

0.880 

1.552 

0.1*1*7 

0.665 

1.71*9 

0.316 

0.521* 

1.858 

0.100 

0.1931 

2.009 

0.0316 

0.0627 

2. Old 

0  to  ±  i 

No  roots 

-1.1831 

9.775 

3.6331 

-1.1251 

9.16 

2.731 

-1.3U21 

6.57 

1.821 

-1.3781 

6.1*2 

1.671 

-1.1*1*91 

6.19 

1.1*21 

-1.5281 

5.97  . 

1.1631 

-1.5811 

5.88- 

1.081 

-1.6311 

U.155 

0.9201 

-1.7321 

l*.03 

0.7651 

-2 .OOOi 

3.78 

1.981 

-2.236i 

3.67 

1.621 

h 

P.  / _ \  2 

2a 

X+2*i 

.61*9 

5.81*3 

.61*7 

5.7U7 

.658 

5.323 

.668 

l*.99l* 

.701 

1*.380 

.732 

3.981 

.761* 

3.726 

.788 

3.635 

.801 

3.62 

.850 

3.662 

.889 

3.80 

.921* 

1*.08 

.968 

l*.l6 

2.27 

5.28 

27.25 

2.71* 

3.73 

28.0 

2.79 

2.65 

19.15 

2.78 

2.51 

19.52 

3.08 

2.22 

20.0 

3.36 

1.81' 

20.1* 

3-UU 

1.75 

20.7 

2.77  • 

1.1*8 

11.0 

3 .01* 

1.17 

10.81 

0.956 

11.76 

10.72 

1.012 

10.5 

10.78 

.11*. 


Because  these  computations  were  performed  in  the  closing  days  of  the 

contract,  there  has  been  scant  opportunity  to  check  them  against  the  data 

of  others.  However,  during  the  course  of  the  contract  the  writer  was  in- 

2 

vited  to  review  a  preliminary  report  of  a  series  of  measurements  on  barium 
titanate  cylinders.  This  report  declares  that  the  lowest  frequency  of 
vibration  in  fully  polarized  cylinders  appears  to  obey  the  following  re¬ 
lationships  in  Which  and  are  constants, 

(■f  D)2  =  M,  [ 1  -  }  0  1  (22) 

and 

(f  L)1  =  N,  [ !  -  »*  (jr)’]  (23) 

where  D  and  I  are  2a  and  h,  respectively,  in  this  writer’s  notation.  Figs. 

U  and  $  are  displays  of  computed  values  arranged  for  qualitative  comparison 
with  Eqs.  (22)  and  (23).  The  latter  of  course  would  be  straight  lines  of 
negative  slope  if  superposed  on  Figs.  I*  and  5. 

There  is  virtually  no  encouragement  derivable  from  this  comparison. 

However,  in  view  of  the  limited  range  of  h/2a  covered  by  the  confutations, 
one  cannot  conclude  for  or  against  the  theoretical  equations,  Eqs.  (U)  and  (5), 
on  the  basis  of  this  comparison. 


^Work  by  G.  K.  Lacey,  Jr.,  Diamond  Ordnance  Fuse  Laboratory,  Washington,  D.  C., 
publication  pending. 
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Section  IT 

DIRECTION  OF  FUTURE  WORK 

Three  areas  within  the  search  for  other  solutions  have  been  designated  as 
incompletely  tested.  These  are  (1)  the  assumption  that  pj  and  Bj  are  complex, 
(2)  the  discovery  of  additional  solutions  to  the  wave  equation,  and  (3)  a  more 
comprehensive  evaluation  of  Eqs.  (U)  and  (5)  for  j  -  0. 

It  is  not  possible  to  say  which  has  the  most  promise.  Further  testing  of 
(1)  could  probably  be  completed  in  the  least  time  and  would  therefore  be  rec¬ 
ommended  for  priority.  It  would  be  the  choice  of  the  writer  to  put  second 
priority  upon  a  part  of  (3),  namely  the  evaluation  of  fundamentals  sufficient 
to  cover  the  full  range  of  h/2a,  before  undertaking  (2).  The  evaluation  of 
overtones  would  be  designated  last,  it  being  contingent  upon  earlier  results. 

The  further  approach  to  complex  propagation  constants  would  be  illustrated 
by  the  following  frequency  and  displacement  equations,  wherein  the  asterisk 


represents  complex  conjugate. 


T -  ff 


:*+  my  =  f>j  +  »j 


^  <*> 

u*  ’  .4JJ‘(rrXA wt 

Uz  -  it  Sin*}*  +  Bjj-H  Sinfn j+iZ] 

im’*1  \  +  •  %  n-  flVih S  <26> 

+  'J»(Pj'i)[DjjS/n  myz  +  £>]/■*■•  Sin  m^+iZjVe 
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It  is  felt  that  there  are  sufficient  unknowns  in  Eqs.  (2l*)-(26)  to  equal 
the  number  of  equations  if  the  arbitrary  constants  C  and  D  are  not  complex 
conjugates  of  A  and  B  respectively. 

The  approach  to  testing  additional  solutions  of  the  wave  equation  for 
efficacy  in  solving. the  cylinder  problem  is  not  completely  planned  at  this 
time.  A  first  step  would  be  to  assemble  three,  say,  pairs  of  displacement 
equations  that  would  be  next  in  increasing  order  of  complexity  beyond  Bqs. 
(1U)-(19).  Since  these  six  pairs  are  independent  solutions,  there  are  ini¬ 
tially  as  many  arbitrary  constants  as  there  are  pairs  chosen.  The  con¬ 
straints  of  the  problem  would  then  be  applied,  hoping  that  the  conditions 
could  be  satisfied  and  the  problem  solved. 

If  not,  a  second  step  is  the  search  for  a  generating  function  by  Which 
the  six  pairs  and  all  possible  pairs  can  be  predicted.  The  sum  of  all  pos¬ 
sible  pairs  would  then  be  tried  as  a  solution. 

If  the  second  step  fails,  a  third  step  is  a  renewed  assault  upon  the 
manipulation  of  the  general  equations  which  aye  Eqs.  (20)  and  (21).  For 
this  there  are  no  definite  guidelines  known  at  this  time. 
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APPENDIX  B 
A  ddi.iion.al  Solutions 
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be  the  radial  and  axial  components  of  shear  d  is  placement , 
respectively ,  loifk.  time  Suppressed.  In  order  of  increasing 

decree  in  A  and  z  ,  the  di's placements  are 
UAJ  -  fAj  J,  cos 

and 

Uzd  =  *nAa  Jo  *'"■  » 

t/„a 

»  -  p  A,  J.  «''«■  ’ 


and 


and 


=  pBj[m(tnn  Jo  cos  +  f>*  Jt  sin)  +  f>J,  cosj 
Uzj  -  -  iw6d[r*  (niA  J,  sin  -+ pz  J9  cm)J  f 

a.*  B,[f  nt  (i«aJo  c«5  -+•  pz  J,  sin.)  -  (p%+2mx)  J,  c«jJ 
Uzs  =  B4  [  pa  (  rnA  J,  Sih  ■+  fz  J.  coa)J  ; 


A  r<^ 


u.mtaa 


.4s  pA.  and  m-z.  omit  ted  for  brevity  , 


B  I 


and 


and 


=  fx  Cj  {  **a  [  (p  z)*+  (hnA)2]  J,  cos  -  m*p  a.  Je  (cos +2  mi  s/'/i) 

-  2  p*  J,  (  c  65  +  mz  sin  )  | 
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L>zs  = 
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